Draft version April 4, 2012 

Preprint typeset using I4TgX style emulateapj v. 5/2/1 1 



ON THE DETECTABILITY OF DUAL JETS FROM BINARY BLACK HOLES 

Philipp Moesta'-*, Daniela Alic', Luciano Rezzolla'-^, Olindo Zanotti', Carlos Palenzuela^-^ 

Draft version April 4. 2012 

ABSTRACT 

We revisit the suggestion that dual jets can be produced during the inspiral and merger of supermassive black 
holes when these are immersed in a force-free plasma threaded by a uniform magnetic field. By performing 
independent calculations of the late inspiral and merger, and by computing the electromagnetic emission in a 
way which is consistent with estimates using the Poynting flux, we show that a dual-jet structure is present but 
energetically subdominant with respect to a non-collimated and predominantly quadrupolar emission, which is 
similar to the one computed when the binary is in electrovacuum. While our findings set serious restrictions 
on the detectability of dual jets from coalescing binaries, they also increase the chances of detecting an EM 
counterpart from these systems. 



1. INTRODUCTION 

The inspiral and merger of supermassive black holes (BHs) 
will represent a secure source for the planned space-borne 
gravitational- wave (GW) detectors. Together with the GW 
sign al, this process is e xpected to be acc ompanie d either be- 
fore'Bode e t"aL|(|20n| ); |Farris et al.| ( |20ri ) or after Rossi etal.' 
pOlO); Zanotti et al.|(|2010)l the merger by the emission of 
electromagnetic (EM) radiation, thus providing a perfect ex- 
ample of multi-messenger astronomy. Should a "simultane- 
ous" detection take place, this would not only help to local- 
ize the GW source and provide its redshift, but also address 
a num ber questi ons in astrophys ics and cosmology Haiman] 
[etaL] C2009); Sesana et al.'lOTll. 

As the merger between two galaxies takes place and the 
BHs get closer, a circumbinary accretion disk is expected to 
form. Because the radiation-reaction timescale over which the 
binary evolves is much longer than the accretion timescale, 
the disk will slowly follow the binary as its orbit shrinks. 
However, as the binary separation becomes ^ 10^ M, where 
M is the mass of the binary, the radiation-reaction timescale 
reduces considerably, becoming smaller than the accretion 
one. When this happens, the disk evolution disconnects from 
the binary, the accretion rate reduces, and the bin ary performs 
its final orbits in an inner regi on poor of gas |Armitage &| 
[Natarajan (2Q02jl; puetaL] ( p003| l. 

This basic picture represents the astrophysical backdrop of 
a simple model which has been used by a number of authors 
to model the EM emission from the BH binary. More specif- 
ically, assuming that the disk is threaded by a coherent and 
large-scale magnetic field which is anchored to the disk, this 
will also permeate the "evacuated" region where the binary 
is rapidly shrinking and provide an effective way of coupling 
the binary's orbital motion to the generation of an EM sig- 
nal. This scenario has b een considered both for spa c etimes 
in ele ctrovacuum (EV) Palenzuela et al.'('2009','2010^;'Mosta' 
|et al.| ([2010) and for spacetimes filled by a tenuous plasma in 
the force-free (FF) approximation (Palenzuela et al.| ( |20lObja[ ); 

' Max-Planck-Institut fiir Gravitationsphysik, Albert-Einstein-Institut, 
Potsdam-Golm, Germany 

^ Department of Physics and Astronomy, Louisiana State University, 
Baton Rouge. LA, USA 

^ Canadian Institute for Theoretical Astrophysics, Toronto, Ontario 
M5S 3H8, Canada 

TAPIR, MC 350-17, California Institute of Technology, Pasadena, 
CA, 91125 USA 



[Neilsen et aLl ( |2010| ). Reference [Mosta et al.| ( |201 0l), in par- 
ticular, has considered a series of spinning BH binaries and 
studied the dependence of the gravitational and EM signals 
with different spin configurations. All in all, it was found that 
EM radiation in the lowest £ = 2, rn = 2 multipole accurately 
reflects the gravitational one. Furthermore, the efficiency of 
the energy emission in EM waves was found to scale quadrat- 
ically with the total spin and to be ^ 13 orders of magni- 
tude smaller than the one in GWs for realistic magnetic fields. 
However, the prospects of detecting an EM counterpar t have 
become larger when it was pointed out in Palenzuela et al. 
(20 10b I that during the inspiral in an FF plasma, a dual-jet 
structur e forms, as a generalization of the Blandford-Znajek 
process |Blandford & Znajek| ( |1977| l to orbiting BHs. Under 
these conditions, the EM energy can accelerate electrons and 
lead to synchrotron radiation. 



We extend here the investigations made in Mosta et al. 
( 2010 1 by considering the evolution in an FF regime and prop- 
erly computing the EM emission as the net result of ingoing 
and outgoing radiation. We confirm in this way the presence 
of a dual-jet structure, but also find that this is energetically 
subdominant with respect to a non-collimated and predomi- 
nantly quadrupolar emission, at least during the late inspiral 
and merger we have considered. Hence, if the simplified sce- 
nario considered here is realized in astrophysical configura- 
tions, it will be difficult to reveal the emission from the dual 
jets, but it will also be easier that an EM counterpart to binary 
BH mergers can be detected. 

2. METHODOLOGY 

As mentioned above, we assume the magnetic field to be 
anchored to the disk, whose inner edge is at a distance of 
^ 10^ A/ and is effectively outside of our computational do- 
main, while the initial binary separation is only D — % M . 
Although the large-scale magnetic field is poloidal, it is set 
to be initially uniform within the computational domain, 
i.e. B' = (0,0, So) with Bq - IQ-'^/M. We note that 
although astrophysically large, the initial magnetic field has 
an EM energy several orders of magnitude smaller than the 
gravitational-field one, so that the EM fields can be treated as 
test-fields and all the results for the luminosity scale quadrat- 
ically with i?o- The electric field i?* is initially zero, bu t 
it rapidly reaches a consistent solution Mosta et al. (20101. 
As BH initial data we consider binaries with equal masses 
but with two different spin configurations; namely, the sq- 



2 



binary (both BHs are nonspinning) and the sg -binary (both 
BHs have spins aligned with the orbit al angular momentum 



with (J/M^)i = ai = 02 = 0.6; see Rezzolla et al.| ( |2008| l; 



[Reisswig et al.| ( [2(309| l for details). These two contigiirations 
allow us to study both the contributions coming from the bi- 
nary's orbital motion, but also those related to the spins of 
the two BHs. The setup for the numerical grids used in the 
simulations consists of nine levels of mesh refinement with a 
fine resolution of Ax/AI = 0.025. The wave-zone grid has 
a resolution of Ax/M ~ 1.6 and extends from r — 24 M 
to r = 180 A/, while the outer (coarsest) grid extends to 
819.2 M. In our implementati on of the FF equations w e solve 
the same set of equations in Palenzuela et al. ( 2010a i but do 
not enforce "by-hand" the FF condition. Rather, we imple- 
ment a damping scheme which drives the solution to satisfy 
the FF condition and avoids i nconsistent EM fields; we will 
present our approach in Alic & others '(201 li 

The calculation of the EM and gravitational radiation gener- 
ated during the inspiral, merger and ringdown is arguably the 
most important aspect of this work as it allows us to establish 
what are the characteristics of the emissions in two channels. 
For the GW sector, we compute the emission via the Ne wman- 
Penrose (N P) curvature scalar ^'4 as detailed in PoUney et al. 
(20071 and |Mosta et al. (2010jl. In practice, we define an or- 
thonormal basis in the three space {r,6,<p), with poles along 
z. Using the normal to the slice as timelike vector t, we con- 
struct the null orthonormal tetrad {I, n, m, m}, with the bar 
indicating a complex conjugate 

l=^{i + r), n=-^(i-f), m=-^(0 + i0), 

(1) 

by means of which we project the Weyl curvature tensor 
Cap^s to obtain vl'4 = Cap^sn'^fh^n^fh^ ■ For the EM sector, 
instead, we use two equivalent approaches to cross-validate 
our measures. The first one uses the NP scalars $0 (for the in- 
going EM radiation) and $2 (for the outgoing EM radiation), 
defined using the same tetrad Teukolsky ( 1973) 
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It is always useful to remark that by construction quantities 
such as ^4,<i>o,<i>2 are well-defined only at very large dis- 
tances from the sources. Any measure of these quantities in 
the strong-field region risks to be incorrect. As an example, 
the EM energy flux does not show the expected 1 /r^ scaling 
when <l>o, <&2 are rnea sured at distances of r ~ 20 M |Palen-| 
zuela et al. ( 2010b|a 1, which is instead reached only for r > 
100 TW. As we will show later, this results in a considerable 
difference in the estimate of the non-coUimated EM emission. 

Since our choice of having a uniform magnetic field 
within the computational domain has a number of drawbacks 
(e.g. nonzero initial $2 , 'i'o), great care has to be taken when 
measuring the EM radiation. Fortunately, we can exploit the 
linearity in the Maxwell equations to distinguish the genuine 
emission induced by the presence of the BH(s) from the back- 
ground one. Following Teukolsky ( 1973 ), we compute the to- 
tal EM luminosity as a surface integral across a 2-sphere at a 
large distance 

= lini / (|$2 - $2,b|' - |$o - *o.Bp) dn , (3) 

r— >oo J 

where (&2.B and $o.b are the values of the background scalars 
induced by the asymptotically uniform magnetic field solu- 



tion in the time dependent spacetime produced by the binary 
BHs. The choice of $2,b and $o,b will be crucial to measure 
correctly the radiative EM emission. 

The background values of the NP scalars $2,67 'I'o,b to be 
used in Equation Q are themselves time dependent and can- 
not be distinguished, at least a priori, from the purely radiative 
contributions. A first prescription takes the background val- 
ues to be time independent and those at the initial time 
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A second prescription involves instead the removal of those 
multipole components of the NP scalars which are not radia- 
tive, namely all those associated with the m = multipoles 
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where {^2,o)e.m=o refer to the to = modes of the multi- 
polar decomposition of $2.0 (^ < 8 is sufficient to capture 
most of the background). Note that because the m = back- 
ground is essentially time independent (after the initial tran- 
sient), the choice (|5]l is effectively equivalent to considering 
as background the final values of the NP scalars as computed 
in an EV calculation of the same binary system. 

While apparently different, expressions (|4]) and (|5]l lead 
to very similar estimates and, more importantly they have a 
simple interpretation in terms of the corresponding measures 
made with. We recall, in fact, that the EM luminosity is 
customarily computed via the integral over a 2-sphere of the 
Poynting flux Si — ^eij^E^ , which is again just the flux 
of the stress-energy tensor as measured now by observers on 
the spatial hypersurface. Of course also this measure is ad- 
equate only far from the binary and it suffers equally from 
a background non-radiative contribution. However, because 
of the linearity in the Maxwell equations, the non-radiative 
contributions can also be removed by introducing background 
values of the EM fields £^g, -Bg and computing the Poynting 

vector as 5".^ = ^€^.j,,{E' - E^){B^ ~ B^). In this con- 
text, then, expressions Q and (|5f correspond respectively to 
setting = 0, B| = = 0) and = El^=o, = 
B^m=a- ^^^s verified that the measures of the EM lu- 
minosity obtained using Equation Q or (|5]l reproduce well 
the corresponding ones obtained with the Poynting flux. As 
a final remark we note that the EM flux in Equation ([3]) is 
not always positive everywhere on the 2-sphere. The negative 
contributions, however, average to zero over one orbit and do 
not represent a radiativ e field. This point was remarked in 
Palenzuela et al. (2010'), where a toy model within the mem- 
brane paradigm was used for the binary. 

3. RESULTS 

Using expression (|3]l with either the prescription (|4|i or (|5]l 
we find that the EM radiation generated during the late inspi- 
ral and merger does contain a dual-jet structure, but also that 
the total energy flux is dominated by a non-collimated emis- 
sion of quadrupolar nature. This is shown in Figure [T] which 
reports snapshots at different times of the EM energy flux on 
a 2-sphere of radius r ~ 100 M for the SQ-binary. The top 
row shows the fluxes as computed with Equation (|4]) and with 
a (non-radiative) component that enhances two lobes of the 
signal and decreases the other two (the first, second and third 
panel refer to the inspiral, merger and ringdown stages, re- 
spectively). Although there are regions with negative values 
which will average to zero over one orbit, the corresponding 



100 
75 
[M] 50 
25 


-100 



1 .Oe-08 
7.8e-09 
5.6e-09 
3.4e-09 
1 .2e-09 
1.0e-09 



-60 



-20 



t=251 M 



594 M 



100 




100 



100 



t=594 M 




100 



100 



100 



condition j4j and via jsj, respectively. Note that 




Fig. 1 . — Snapshots of the EM energy flux during the inspiral and merger of the binary sg. The top and bottom rows show the fluxes as measured via j3j with 

sg as computed through Equation (|5]). Both the coUimated 
and the non-colUmated emission are very similar, ahhough 
there is a 50% enhancement of the radiation in the spinning 
case, both in the total and in the collimated emission. This is 
because most of the radiation is produced by the interaction 
between the BH orbital motion and the background magnetic 
field. Indeed we find the emission in the EV evolution to be 
compar able to the FF one (th is is different from what reported 
in Pale nzuela et al.| ( |2010b|al )). 

A more quantitative assessment of the different contribu- 
tions is shown in Figure [3] which reports the evolution of the 
EM luminosity of the binary sq (left panel) and of the spinning 
binary sg (right panel). Shown with (blue) solid lines is the 
total luminosity as computed through Equation (|5]), while the 
(red) dashed lines refer to the luminosity integrated over a po- 
lar cap with a half-opening angle of 5° and thus representative 
of the emission from the two jets. Also shown with (magenta) 
dotted lines are the luminosities as computed through Equa- 
tion the differences are small and hardly visible for the 
collimated part. 

An accurate measure of the evolution of the collimated/non- 
collimated components is crucial to predict the properties of 
the system when the two BHs are widely separated. This, 
in turn, requires a reliable disentanglement of the collimated 
emission from the non-collimated one and from the back- 
ground. The non-collimated emission measured with Equa- 
tion (|5]l matches well the growth expected if the EM emis- 
sion is mostly quadrupolar and hence with a dependence that 
is the same as the G W one and, at the lowest order, scales 
as L!!?"-™" w ^^"/'^ [Palenzuela et al.|(|2010|i;|M5sta etaL 
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Fig. 2. — Map views of the EM energy flux computed using |3j with j4j. 
The top row refers to the binary so, while the bottom one to the binary sg 
and has a larger emission from the dual jets. 

EM fields could nevertheless induce motion in the plasma and 
secondary radiation which we cannot account for here. The 
second row refers to fluxes as computed with Equation (|5]l, 
where the time dependent background solution is almost com- 
pletely subtracted, so that only the (positive) radiative part re- 
mains, showing a fairly symmetric four-lobe structure. Within 
both approaches, the extended quadrupolar distribution is ac- 
companied by the presence of dual jets during the inspiral, 
and of a single jet from the spinning merged BH. The en- 
ergy flux in the jets is essentially the same in the two cases 
(both during the inspiral and after the merger), but that in the 
non-collimated part is different. In spite of this, the total lu- 
minosities are very similar, as most of the differences cancel 
when integrated over one orbit (cf. Figure |3]l. However, be- 
cause the Poynting-flux structure is different in the two mea- 
sures (Equation Q has a non-radiative part missing in Equa- 
tion (|5|), it could lead to a different secondary emission as the 
EM fields interact with the plasma; this emission cannot be 
investigated within an FF approach but deserves attention. 

The corrections introduced by the spin of the BHs are 
shown in Figure l2] which reports map views (i.e. projections 
on the {x, y) plane) of the EM emission for the binaries sq and 



(20101. On the other hana, a smaller scaling is founa wnen 
the approach Equation (|4]) is adopted; at present it is difficult 
to determine which one is the most reliable scaling. At the 
same time, the frequency evolution for the collimated emis- 
sion coincides in the two approache s, but it does not follow the 
scaling L'^°^^ 



(2010b I 

6/^ 



suggested in Palenzuela et al. 
for boosted BHs; rather, it scales as L™'^ < rj^/^-*^/^. We 

' ' EM ^ 

suspect the accelerated motion of the BHs to be behind this 
difference and longer simulations are needed to draw robust 
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Fig. 3. — Evolution of the EM luminosity at 100 M for the binary sq 
(left) and for the spinning binary sg (right), when AI = 10* Mq and 
-Bo = 10* G. Using jsj, (blue) solid lines show the total luminosity, while 
(red) dashed lines refer to the luminosity in anolar cap of 5°. Shown with 
(magenta) dotted lines are the measures with E); note the presence of a small 
eccentricity for the binary Sfy. 

conclusions. 

Impact on detectability. Assessing the detectability of the EM 
emission discussed abov e is of course of great importance 
and using th e results of Palenzuela et al.| (|2010b|all, IKaplanI 
|et al.| ( pOlT] l have estimated that the short timescales associ- 
ated with the merger will limit the detectability in the radio 
band to less than lyr^^. As reported in Figure [5] the peak 
of the collimated emission is ^ 100 times smaller than that 
of the total emission, making the detection of the dual jets at 
the merger unlikely. Unfortunately, even if the energy flux 
is ~ 8 — 2 times larger near the jets, the lack of knowledge 
about the Lorentz factor of the reprocessed plasma does not 
allow us to say whether the beaming in the jet will be larger 
than that in the extended emission and thus help its detection. 
That said, because the total luminosity at merger is ^ 100 



times larger than in Palenzuela et al. ( 2010b i (mainly because 
we measure it at large distances where it has approximately 
reached its asymptotic value), the detection shoul d be overall 
more likely if most of the assumptions in Kaplan et al.| ( |20TT] l 
are verified. In addition, the dual jets emission could be dom- 
inant in the early inspiral (especially if the BHs are spinning). 
Assuming the scalings for the collimated and non-collimated 
emissions are different and go respectively like ^^^/^ and 
^-^6/3 (^Qji/^Y the two components of the luminosity would 
become of the same order at a separation of ^ 16 (24) M, 
but obviously smaller (the collimated will be smaller by a 
factor ~ 10). Determining more precisely when and if this 
happens requires an accurate frequency scaling which is not 
available yet. As a note of caution we stress that luminosities 
-^EM ^ 10^^ '^''S ^l^o typical of radio-loud galaxies, 

and the determination of an EM counterpart can be challeng- 
ing if such sources are near the candidate event. 

4. CONCLUSIONS 

We have investigated the suggestion that dual jets can be 
produced during the inspiral and merger of supermassive BHs 
immersed in an FF plasma threaded by a uniform magnetic 
field. We have found that the energy flux does have a dual- 
jet structure but is predominantly quadrupolar, with the non- 
collimated emission being about 10—100 times larger than the 
collimated one. Our findings set restrictions on the detectabil- 
ity of dual jets from coalescing BH binaries, but also increase 
the chances of detecting an EM counterpart for astrophysical 
conditions similar to those in this simplified scenario. 
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